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SECTION  I 

INTRODUCTION  AND  SUMMARY 
1.1  BACKGROUND  AND  OBJECTIVES 

The  survivablil ity  of  any  candidate  M-X  basing  system  must 
be  evaluated  using  estimates  of  the  hostile  environment  produced 
by  nuclear  bursts.  Experimental  and/or  numerical  simulations 
can  be  used  to  predict  the  physical  characteristics  of  the 
hostile  environment.  These  predictions  and  their  associated 
uncertainties  can  then  be  used  in  analysis  of  the  structural 
response  of  important  elements  in  the  basing  system. 

This  report  is  concerned  with  one  of  the  physical  processes 
which  affects  the  design  environment  in  the  buried-trench 
multiple  aim  point  (MAP)  basing  concept1:  After  a  nuclear  burst 
within  the  trench  or  on  the  surface  above  a  trench,  a  strong  shock 
will  be  driven  down  the  trench.  The  high  velocity  flow  of  shocked 
air  in  the  trench  will  vaporize,  melt,  crush,  and  mix  trench  wall 
material  (e.g.,  concrete,  soil).  Much  of  the  debris*  will  inter¬ 
act  and  become  entrained  into  the  flow,  and  this  debris  will 
attenuate  the  air  flow  velocity  and  shock  strength  by  absorbing 
momentum  and  energy  through  drag  and  heat  exchange.  Our  objective 
has  been  to  examine  the  degree  of  this  attenuation  and  the  effects 
of  debris  entrainment  ( and  associated  uncertainties )  on  the 
pressure  and  impulse**  delivered  to  the  plug  in  an  MX  trench. 


We  will  use  the  term  "debris”  to  include  both  vapor  and 
condensed  phase  trench  materials  which  become  entrained 
into  the  high  temperature  air  flow. 

Impulse  is  delivered  to  the  plug  both  by  air  pressure  and  by 
debris  impact,  i.e., 

r  2Mflu, 

I  =  /  Pdt  +  a-  d 

J  A 

o 
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1 . 2  APPROACH 

The  1^-D  version  of  the  multi-phase  DICE  code*  was  used  to 
analyze  the  flow  of  hot  air  and  debris  in  covered  MX  trenches, 
and  to  examine  the  effects  of  debris  entrainment  upon  the  loading 
imposed  on  a  rigid  plug  in  the  trench  1800  feet  from  a  1  Mt  on¬ 
line  surface  burst.  For  all  the  cases  which  were  analyzed,  it 
was  assumed  that  the  burst  coupled  6.2  kt  (all  internal  energy) 
to  the  air  in  the  trench  in  the  first  18  feet  in  both  directions 
from  the  burst  (i.e.,  3.1  kt  on  each  side).  A  14  foot  diameter 
trench  buried  6  feet  in  the  ground  is  the  assumed  configuration. 

These  1^-D  numerical  simulations  treated  or  took  into  con¬ 
sideration  the  phenomena  which  are  shown  in  Figure  1,  either 
through  explicit  calculations  of  the  physical  processes,  or 
through  prescription  by  separate  models. 

The  processes  by  which  material  is  removed  from  the  trench 
wall  and  becomes  entrained  in  the  flow  in  the  trench  involve 
complex  stress  wave,  thermodynamic,  and  turbulent  boundary  layer 
and  mixing  phenomena  occurring  under  conditions  which  have  not 
been  accessible  to  experimental  observation  or  investigation. 
Consequently  there  are  major  uncertainties,  for  any  trench 
materials  or  construction,  in  the  mass  flux  of  debris  off  the 
walls,  in  the  physical  state  of  the  debris  (including  size 
distribution  of  condensed-phase  material) ,  and  in  the  timing  of 
effective  entrainment  (mixing)  of  debris  into  the  trench  flow. 


The  basic  DICE  code2,  is  a  2-D  (axisymmetric)  method  to  model 
the  dynamics  of  formation  of  dust  clouds  following  nuclear 
surface  bursts.  It  is  an  implicit,  compressible,  Eulerian 
finite  difference  code  which  can  treat  the  separate  flow  of 
different  phases  (condensed  and  vapor),  material  species, 
and  particle  sizes.  The  different  groups  and  the  basic 
medium  (usually  air)  are  mutually  coupled  by  aerodynamic 
(drag)  and  thermodynamic  (heat  exchange)  processes. 
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In  the  current  program,  we  examined  the  effects  of  these 
uncertainties  on  pressure  attenuation,  trench  wall  recession, 
and  the  plug  environment.  This  was  done  by  numerical  simulations 
in  which  several  substantially  different  models  were  prescribed 
for  the  mass  flux  of  debris,  its  size  distribution,  and  the  timing 
of  its  entry  into  the  trench  flow. 

1.3  SUMMARY  AND  CONCLUSION 

A  baseline  solution  of  the  trench  air  flow  with  no  debris 
(air-only)  produced  a  maximum  pressure  of  6800  bars  on  the  plug, 
and  an  impulse  of  24  bar-sec.  Eleven  calculations  were  made  of 
the  trench  flow  with  various  debris  entrainment  models.  In  all 
cases,  the  maximum  pressure,  and  the  impulse  on  the  plug  (which 
includes  the  contributions  of  both  pressure  loading  and  debris 
mass)  were  substantially  lower  than  in  the  air-only  case.  However, 
the  magnitude  of  the  reductions  was  strongly  affected  by  the  en¬ 
trainment  parameters. 

Three  aspects  of  debris  entrainment  were  examined: 

o  Entrainment  delay  (mixing  length) .  Wall  removal  at 
a  given  station  in  the  trench  begins  as  soon 
as  the  shock  passes.  However,  some  interval  of  time 
(and  therefore  of  distance  behind  the  shock)  elapses 
before  debris  is  mixed  with  the  flow  over  the  entire 
cross-sectional  area  of  the  trench.  To  take  this 
delay  into  consideration,  we  used  the  entrainment 
models  in  the  following  sketch  to  prescribe  the  behavior 
of  the  wall  recession  thickness,  T,  as  a  function  of  local 
gas  flow  velocity,  Ug,  and  the  distance  behind  the  shock 
front,  X. 
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Delayed-Step  Entrainment 
T  =  g(x)  is  the  wall 


recession  rate 


Ramp  Entrainment 
from  debris  entrainment 


The  entrainment  rate  parameter,  B,  is  the  maximum  ratio 
of  wall  recession  rate,  T ,  to  the  local  velocity  of  the 
gas  flow,  Ug. 

Plug  loading  was  found  to  be  very  sensitive  to  the  en¬ 
trainment  delay.  Where  significant  entrainment  occurred 
near  the  shock  (say  20-50  feet  behind  the  shock),  plug 
loading  was  reduced  by  orders  of  magnitude.  For  example 
entrainment  of  the  debris  from  0.15-in.  of  wall  recession 
(with  entrainment  starting  18  feet  behind  the  shock) 
reduced  the  maximum  pressure  on  the  plug  from  6800  bars 
to  78  bars  and  the  impulse  from  24  to  3  bar-sec  (as  com¬ 
pared  to  the  air-only  case)  . 

o  Entrainment  rate  (and  entrained  mass) .  Increasing  the 
entrainment  rate  (and  mass)  reduces  loading  on  the  plug; 
it  may  be  possible  to  effectively  quench  the  gas  flow 
through  techniques  which  enhance  entrainment  within  say 
20-50  feet  behind  the  shook. 
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Figure  2  shows  the  effects  of  entrainment  delay 
(mixing  length)  and  the  entrainment  rate  function, 

3  =  T/U  ,  on  the  pressure- time  history  at  the  plug 

o 

for  several  cases.  The  important  effects  of  an  order 
of  magnitude  increase  in  the  entrainment  rate  parameter 
can  be  seen  by  comparing  Cases  3  and  5.  Similarly,  the 
effects  of  reducing  the  entrainment  distance  by  a  factor  of 
two  (from  108  to  54  feet)  can  be  appreciated  by  comparing 
Case  5  to  Case  6 . 

o  Debris  size.  Analysis  in  which  the  debris  mass  is  initially 
divided  evenly  between  0.1,  1,  and  10  cm  particles  were 
compared  with  cases  using  just  .01  cm  (=  100y)  particles. 
The  smaller  distribution  reduced  plug  peak  pressure  loading 
by  about  a  factor  of  5  more  than  the  larger  distribution. 
Both  drag  and  heat  transfer  depend  on  particle  size;  to 
determine  which  is  more  important,  some  analyses  were  re¬ 
peated  with  no  heat  transfer  being  allowed  between  the  gas 
and  the  particles.  Plug  loading  was  essentially  the  same, 
with  or  without  heat  transfer .  This  indicates  that  the  ad¬ 
dition  of  debris  mass  affects  the  mixed-phase  flow  dynamics 
primarily  through  drag  interactions  rather  than  through 
thermal  interactions. 
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Over  Pressure  on  the  Plug  (bars) 


SECTION  II 
ENTRAINMENT  MODEL 

In  the  multi-phase  DICE  analyses,  passage  of  the  air  shock  is 
assumed  to  initiate  removal  of  wall  material  and  to  start  its  mixing 

with  the  air  flow.  Debris  is  assumed  to  be  solid-phase  silica 
at  the  time  of  mixing,  but  thermal  interactions  can  subsequently 
melt  and  vaporize  the  particles.  The  initial  debris  particle 
size  is  specified  (sizes  from  lOOp  to  10  cm  were  considered) . 

The  entrainment  delay  and  mass  rate  are  specified  by  the 
following  models: 


T  =  wall  recession 
x  =  distance  behind  shock 
Pp  =  2  g/cm^  =  trench  wall  material  density 

$  =  mass  entrainment  rate  p er  unit  area  at  any 
trench  station 


TpP  =  6pPUg 


where  6  = 


=  entrainment  rate  parameter  which  is  a  function  of  x 


U  =  local  gas  flow  velocity 
© 


Note  that  both  3  and  Ug  are  functions  of  the  instantaneous 
distance  of  the  trench  station  behind  the  shock  front. 

The  variation  of  3  with  x  was  described  by  delayed  step* 
function  models  and  by  ramp- function  models.  The  parameters 
for  several  of  the  models  used  are  shown  in  Figure  3,  page  15. 

The  variation  of  Ug  with  x  depends  on  the  flow  dynamics 
behind  the  shock.  These  flow  dynamics  are  affected  by  the 
entrainment  of  debris;  hence  Ug  is  coupled  to  the  entrainment 
delay  and  the  entrainment  mass  rate,  M.  This  coupling  is  taken 
into  consideration  explicitly  in  the  numerical  calculations. 

Models  used  in  the  DICE  calculations  for  other  processes  , 
including  aerodynamic  drag,  thermal  interactions,  and  simulation 
of  tunnel  wall  expansion  and  venting,  are  described  in  the 
Appendix. 
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SECTION  III 
PLUG  LOADING 

3.1  EFFECTS  OF  ENTRAINMENT  DELAY  AND  MASS  LOADING 

Table  1  summarizes  the  time  of  arrival,  maximum  pressure, 

(P^)  and  impulse  (1^)  delivered  to  the  plug  for  six  variations 
in  wall  debris  entrainment  characteristics.  Case  1,  the  limiting 
case  with  no  debris  (air  only),  was  run  for  comparison*.  Cases  2 
through  6  use  the  entrainment  rate  parameters  6 (x),  shown  on 
Figure  3.  Figure  4  lists  the  maximum  pressures  and  impulses 
directly  on  the  6(x)  profiles,  in  order  to  illustrate  the 
sensitivity  of  PM  and  IM  to  uncertainties  in  debris  entrainment. 

A  comparison  of  the  two  delayed  step  -  function  cases  (2  and  4) 
shows  that  introducing  mass  near  the  shock  front  is  very  effective 
in  attenuating  the  plug  P^  and  1^.  When  the  "mixing  length"  (L  in 
Table  1)  was  reduced  from  108  feet  to  18  feet,  the  peak  pressure  on 
the  plug  dropped  from  2500  bars  to  78  bars  (factor  of  32)  and  the 
impulse  dropped  from  12  bar-sec  to  3  bar-sec  (factor  of  4).  Simi¬ 
larly,  results  when  ramp  functions  of  8  were  used  (Cases  3,5,  and 
6)  show  a  strong  dependence  on  the  debris  entrainment  character¬ 
istics  near  the  shock  front. 

Figure  5  is  a  plot  of  impulse  and  time  of  arrival  versus 
maximum  pressure  on  the  plug  for  Cases  1  through  6.  The  impulse 
delivered  to  the  plug  varies  roughly  with  the  .6  power  of  maximum 
pressure,  IM  ~(PM)-6,  and  the  time  of  arrival  varies  roughly  with 
the  -.3  power,  TOA  ~(P^)"’3.  The  generality  of  these  relationships 
is  not  known.  However,  if  validated,  these  are  significant 


Results  from  the  MARVEL  test  (2.1  kt  burst  in  a  3.2-foot  dia. 
tunnel)  have  shown  that  the  air-only  model  is  a  poor  represent¬ 
ation  of  what  actually  happens.  Some  incorporation  of  wall 
ablation  and  scouring  is  needed  to  explain  the  MARVEL  observ¬ 
ations  . 
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Table  1. 


Maximum  Pressures  and  Impulses  Delivered  to  Plug  for  Variations 
in  Debris  Entrainment  Model. 


In  all  these  debris  models,  total  mass  was  Initially  divided  equally  between 
0.1,  1,  and  10  cm  particles. 

IS 


TT 


Distance  Behind  Shock  Front  (ft)  —  x 


Figure  3.  Mass  Entrainment  Parameter,  6,  versus  Distance  Behind 
Shock  Front,  x. 
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Figure  4 .  Maximum  Pressure  and  Impulse  Delivered  to  the  Plug  for  Variations 
in  Mass  Entrainment  Parameter,  0  *  B(x). 


relationships  in  that  they  indicate  debris  uncertainties  which  lead 
to  factors  of  10  variations  in  Pm  produce  only  factors  of  -4  un- 
certainty  in  1^. 

The  time-resolved  pressures  on  the  plug  and  the  cumulative 
impulse  delivered  to  the  plug  are  shown  in  Figure  6  for  two  cases. 
Cases  2  and  5  differ  in  the  treatment  of  the  108  feet  behind  the 
shock  front.  When  a  ramp  dependence  for  8(x)  (Case  5)  is  used  in 
place  of  a  delayed  step  function  (Case  2)  the  following  related 
effects  are  observed: 

o  the  maximum  pressure  is  reduced  by  a  factor 
of  -150. 

o  the  total  impulse  delivered  to  the  plug  is 
reduced  by  a  factor  of  -15. 

o  the  time  of  arrival  increases  by  a  factor  of 
-4  to  5. 

o  the  duration  of  the  pressure  pulse  on  the  plug 
increases  by  a  factor  of  -4  to  5. 

Figure  2  summarizes  the  sensitivity  of  the  pressure  loading 
to  variations  in  the  wall  debris  entrainment  rate  parameter  6(x). 
The  pressure- time  histories  have  been  shifted  to  begin  at  the  times 
of  arrival  in  Table  1.  The  important  effects  (-  factor  of  40  in 
P^)  of  an  order  of  magnitude  change  in  B=T/Ug  can  be  seen  by 
comparing  Case  3  ("6  =  10"^)  and  Case  5  (s’  =  10'^).  Similarly, 
the  effects  (-  factor  of  4)  of  reducing  the  mixing  length  L  by  a 
factor  of  2  can  be  appreciated  by  comparing  Case  5  (L=108  feet) 
to  Case  6  (L=54  feet) .  The  plug  loading  characteristics  (P  vs 
time,  PM,  and  IM)  are  all  sensitive  to  the  entrainment  of  wall 
debris  material  near  the  shock  front. 
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Figure  6.  A  Comparison  of  Pressure  and  Impulse  Time  Histories  using  Two 

Wall  Debris  Entrainment  Models.  (Plug  is  Located  1800  ft  from  a 
1  Mton  On-Line  Surface  Burst.) 
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Impulse  (bar-sec) 


3.2 


EFFECTS  OF  PARTICLE  SIZE 


The  multi-phase  mixture  of  air,  solid-liquid  debris  particles, 
and  debris  vapor  interacts  through  drag  and  thermal  mechanisms. 

The  drag  forces  act  on  both  the  air  and  debris;  these  forces 
generally  accelerate  the  debris  material  while  decelerating  the 
high  velocity  air.  The  thermal  interactions  involve  exchanges 
of  energy  primarily  through  ablation.  Both  drag  and  thermal 
interactions  are  a  strong  function  of  solid/liquid  particle  size. 


The  influence  of  uncertainties  in  the  particle  size  distri¬ 
bution  is  shown  on  Figure  7.  The  only  difference  between  the  two 
numerical  simulations  is  that  Case  7  uses  a  single  debris  particle 
size  of  D=100vi  =  .01  cm  as  compared  to  the  nominal  particle  size 
distribution  of  D=.l,  l,  and  10  cm  particles*  used  in  Case  3. 


Figure  7  shows  that  if  wall  debris  is  assumed  to  mix  with 
the  shocked  air  in  "small"  particles  (D  =  100u  =  .01  cm),  as 
compared  to  the  nominal  particle  size  distribution  (D  =  .1,  1, 
and  10  cm) ,  then  the  enhanced  drag  and  thermal  interactions  will 
cause  the  maximum  pressure  to  drop  by  an  additional  factor  of  -5 
and  the  impulse  to  drop  by  an  additional  factor  of  -2. 


In  order  to  evaluate  the  relative  importance  of  drag  as 
compared  to  thermal  interactions,  two  numerical  simulations 
(Cases  8  and  9)  were  performed  which  ignored  all  thermal  mech¬ 
anisms,  i.e.  no  ablation  or  conduction  was  allowed  (similar  to 
Cases  3  and  7  on  Figure  7  except  there  are  no  thermal  inter¬ 
actions).  Table  2  summarizes  the  results  of  these  numerical 
simulations . 


*  In  the  nominal  case,  the  mass  is  initially  divided  equally 
among  each  particle  size  group,  i.e  1/3  of  the  mass  is  in 
.1  cm  particles,  1/3  in  1  cm  particles,  and  1/3  in  10  cm 
particles.  This  nominal  particle  size  distribution  is  used 
in  all  cases  except  Cases  7  and  9. 
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ressure  at 


Impulse  at 


TABLE  2.  Maximum  Pressures  and  Impulses  Delivered 
to  the  Plug  for  Variations  in  Particle 
Size  Distribution  and  Thermal  Interaction 


Case 

Thermal  Model 
Active? 

Particle 

Size 

Distribution 

DP 

(cm) 

Maximum 

Pressure 

Pm 

Max imum 
Impul se 

h\ 

3 

Yes 

.1,  1,  and  10* 

621 

6.4 

7 

Yes 

.  01 

107 

2.4 

8 

No 

.  1 ,  1 ,  and  10  * 

690 

6.9 

9 

No 

.01 

136 

3.7 

*  Mass  is  divided  equally  among  .1,  1,  and  10  cm  particles 
in  these  cases. 


Tcl ble  2  indicates  that  the  in  -  trench  mixed  phase  flow 
dynamics  are  dominated  by  araj  interactions  as  opposed  to 
thermal  interactions.  The  maximum  pressure  is  changed  by  a 
factor  of  only  -1.2  when  the  thermal  model  is  turned  off  and 
the  initial  particle  size  distribution  is  held  constant  (compare 
Case  3  to  Case  8  and  Case  7  to  Case  9).  Thus,  the  entrained 
wall  debris  mass  and  particle  size  distributions  are  important 
in  determining  the  plug  loading  characteristics  (PM,  IM) . 


SECTION  IV 

IN- TRENCH  RESPONSE 


I 


z=o 


///////////////////// / m  mn 


Shock  Pressure,' 
P 


at  Range  Z 


Soil 

Overburden 


Trench 


Tf™  rp/rrrrnm illinium  in / /_/  n inn  n urn 

z  - — *•  z=z 


T  =  wall  recession  caused  by  Z=  In-trench  distance  from  burst 

debris  becoming  entrained 
in  the  air  flow 


4.1  SHOCK  PRESSURE  ATTENUATION 

The  shock  pressure  attenuation  curves  are  obtained  from  the 
numerical  simulations  by  plotting  the  shock  pressure  and  assoc¬ 
iated  range  at  a  series  cf  times,  as  sketch  below: 
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Peak  and  Shock 
pressure 

3=P(Z),  shock  pressure  attenuation  curve 

•Pressure-distance  profiles  at 
/arious  times, 


Z  =  In-trench  distance  from  burst 


Note  that  the  peak  pressure  and 

the  shock  pressure  are 
not  necessarily 
identical,  e.g.  see 
Fig.  8 
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Figure  8  shows  the  shock  pressure  attenuation  curves  and 
examples  of  predicted  pressure  profiles  for  two  debris  entrain¬ 
ment  lengths  of  L  =  18  feet  and  108  feet  using  a  step  entrainment 
function  with  6  =  10‘3  (Case  10  vs  Case  2).  Note  that  the 
attenuation  in  the  two  cases  is  very  different.  Also,  the 
qualitative  nature  of  the  pressure  profiles  is  different.  When 
L  =  18  feet,  the  additional  entrained  mass  near  the  shock  front 

causes  the  peak  pressure  to  (1)  decay  much  more  rapidly  as  compared 
to  the  L  -  108  feet  case,  and  (2)  occur  behind  the  shock  front. 

Figure  9  shows  the  effect  of  varying  the  mass  loading 
rate  parameter,  8,  for  the  step  entrainment  function  with  L  =  18 
feet.  This  figure  shows  the  shock  pressure  attenuation  curves 
for  8  =  .01,  .001  ,  and  .0001. 

Figure  10  indicates  the  effect  of  wall  debris  particle  size 
on  the  shock  pressure  attenuation  characteristics.  In  most  of 
the  cases  performed  in  this  study,  the  wall  debris  was  assumed 
to  enter  the  air  flow  with  1  mm,  1  cm,  and  10  cm  particle 
diameters.  The  three  sizes  were  evenly  distributed  according  to 
mass;  however,  on  Figure  10,  two  cases  are  shown  where  the  wall 
debris  was  assumed  to  consist  of  only  100  micron  (.01  cm)  parti¬ 
cles.  As  Figure  10  shows,  the  100  micron  particles  are  very 
effective  in  attenuating  the  shock  pressures;  this  attenuation 
is  primarily  due  to  particle-air  drag  interactions.  Thus,  "fast" 
mixing  of  small  particles  with  the  shock  air  will  be  effective  in 
attenuating  the  shock  pressures  in  the  trench. 

4.2  WALL  RECESSION 

Figure  11  shows  two  comparisons  of  wall  recession  thick¬ 
ness  versus  distance  profiles.  In  the  upper  portion  of  this 
figure,  (Figure  11a)  the  wall  recession  versus  distance  profiles 
are  compared  for  8  =  10*^  and  F  =  10"^  with  a  step  entrainment 
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Figure  9.  Shock  Pressure  vs  Distance  for  Various  Wall  Debris 
Entrainment  Parameters,  F*.01,  .001,  and  .0001 
(L=18  ft,  8  is  a  step  function). 
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Figure  11a.  Wall  Recession  vs  Distance  Profiles  for  g  =  10  ^ 

and  "B  =  10"^  (L=18  feeL  with  Step  Loading  Function). 
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Figure  lib.  Wall  Recession  vs  Distance  Profiles  with  Entrained 
Particle  Distributions  of  Dp  ■  .01  cm  and  D  «*.l,  1, 

10  cm  (F-10"4  and  L-108  ft  with  Linear  Loading  Function) 
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Wall  Recession  (cm) 


function  and  L=18  feet.  Note  that  even  the  fraction  of  an  inch 
recession  (T-.15  in.)  in  the  6  =  10'4  case  (Case  4)  causes  the 
maximum  pressure  at  the  plug  to  drop  to  80  bars  as  compared  to 
the  -6800  bars  in  the  air-only  case  (Case  1). 

Figure  lib  compares  the  wall  recession  profiles  for  a 
variation  in  debris  particle  size  (Case  7  vs  3 ) .  In  Case  7, 
the  wall  debris  particle  diameters  are  entirely  100  microns  = 
.01  cm.  In  Case  3,  the  nominal  particle  size  distribution  of 
Dp  =  .1,  1,  and  10  cm  is  used  for  the  entrained  wall  debris. 
Note  that  in  the  100  micron  case,  less  mass  has  been  scoured 
off  the  walls,  but  the  attenuation  with  distance  is  greater, 
as  previously  shown  on  Figure  10. 
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APPENDIX  A 


PHYSICAL  MODELS  USED  IN  1%-D  DICE  CALCULATIONS 

A  1%-D  calculation  can  be  used  to  simulate  a  flow  problem 
that  satisfies  two  conditions: 

1.  There  is  a  "well-defined"  principal  flow 
direction . 

2.  Transfer  of  mass,  momentum,  and  energy  trans¬ 
verse  to  the  principal  flow  direction  can  be 
modeled  so  that  numerical  integrations  are 
only  performed  in  the  principal  flow  direction. 

Within  the  context  of  the  present  parametric  analyses  of 
mixed  phase  flow  in  trenches,  these  conditions  are  adequately 
satisfied  during  most  of  the  in-trench  flow  dynamics.  These 
conditions  are  not  satisfied,  however,  when  the  trench  roof 
displacement  approaches  the  trench  radius  in  size.  Therefore, 
the  1%-D  analyses  are  probably  poor  simulations  near  the  burst 
point. 

The  physical  models  described  in  this  section  are  of  two 
types : 

1.  General  Physical  Models 

Drag  Model 

Heat  Exchange  Model 

2.  Models  for  Treating  Effects  of  Transverse  Flow 
Interactions . 

Venting  Model 

Wall  Expansion  Model 

(The  debris  entrainment  model,  which  falls  in  this  group, 
is  described  in  Section  II  of  the  text.) 
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A .  1  AERODYNAMIC  DRAG 


In  DICE  mixed  phase  flow  calculations,  the  dirt  particles 
and  air  experience  mutual  drag  interactions.  The  basic  model 
used  for  the  DICE  drag  interaction  force  (Fp)  on  a  particle  of 
diameter  (Dp)  moving  with  velocity  (Vp)  in  a  gas  stream  of 
density  (Pg)  and  velocity  ( V g )  is: 

2 

^P_ 

4 


-*  irDr 

n  ’  ■*»,  CVV1  vp-vg  I  cDC — 1—3 


(A.l) 


where  the  drag  coefficient  (Cp)  is  defined  by: 

CD  =  0<6  +  H 


(A. 2) 


and  Re  is  the  particle  Reynold's  number, 


Re 


_  V  VVg  IPg 


(A. 3) 


where 


U  =  (2.42T)0’78  x  10"6  poise,  T  in  °K 


The  relative  masses  of  particles  and  gas  in  a  computational 
cell  are  used  to  guarantee  that  changes  in  the  dirt  particle 
momentum  are  compensated  by  oppositely-directed  changes  in  the 
air  momentum  (see  Reference  2  for  additional  details) . 

A. 2  THERMAL  INTERACTION 

The  DICE  thermal  interaction  model  couples  dirt  particle 
ablation,  condensation,  and  internal  heat  conduction  with  changes 
in  the  air  internal  energy.  Energy  is  exchanged  between  the  gas 
phase  (air  and  dirt  vapor)  and  condensed  dirt  particles  by  black 
body  radiation  and  convection.  The  convective  exchange  model 
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includes  the  effects  of  the  strong  radial  gas  flow  that  accom¬ 
panies  small  particle  ablation  or  condensation.  for  a  complete 
discussion  of  the  DICK  thermal  interaction  model,  refer  to 
Appendix  C  of  Reference  4. 

Sample  radiative  and  convective  extinction  times  for  dirt 
particles  at  a  gas  temperature  of  20,000°K  (also  assuming 
P=1  kbar,  Re=104,  and  a  small  dirt  to  air  mass  ratio)  are  shown 
in  Figure  12.  An  extinction  time  is  the  time  required  for  a 
dirt  particle  of  a  given  initial  diameter  to  ablate  completely 
at  constant  gas  temperature.  Observe  that  the  convective  ex¬ 
tinction  time  varies  as  the  square  of  the  particle  diameter 
while  the  radiative  extinction  time  varies  as  the  particle 
diameter.  Thus  the  lOOp  particles  have  radiative  extinction 
times  of  about  300  usee  as  compared  to  3  msec  for  1  mm  particles. 

A.  3  1J*-D  VENTING  AND  WALL  EXPANSION 

Roof  venting  and  trench  wall  expansion  are  coupled  in  the 
lis-D  DICE  solutions.  An  equivalent  square  trench  is  defined 
as  shown  in  Figure  13,  parts  a  and  b.  The  width  (W)  of  a  square 
trench  equivalent  in  cross-sectional  area  to  a  cylindrical  trench 
of  radius  (R0)  is  defined  by 


W  =  /T  R0  (A. 4) 

H  is  the  combined  roof  and  overburden  thickness  of  mean 
density  (ps) . 
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The  models  assume  that  the  roof  displaces  upward  (with 
rigid  body  displacement  Xr)  as  a  trapezoidal  wedge  and  that 
the  trench  gases  vent  through  the  resulting  channel  of  width 
(w) .  The  moving  roof  obeys  the  simple  equation  of  motion. 


(A.  5) 


If  p  is  constant  during  an  integration  cycle,  the  roof 
velocity  and  displacement  are: 


(P'Po) 


t  +  xT 


( A.  6) 


(P'Po)  2 

xR  =  -  t4,  +  xR  t  +  xR 

R  2psH  R°  ° 


( A.  7) 


where  xrq  and  xrq  are  the  roof  velocity  and  displacement  at 
the  beginning  of  the  integration  cycle. 

The  trench  is  modeled  as  completely  open  when  the  vent 
channel  area  equals  the  trench  roof  area  or  when 

Xr  =  W/2cos0  (A. 8) 

Prompt  venting  is  assumed  to  occur  for  all  non-zero 
values  of  xr. 

The  vent  fraction  (n)  is  defined  as  the  ratio  of  vent 
area  to  equivalent  square  trench  roof  area, 


D  -  2xrCos9/W 


(A. 9) 


The  venting 
defined  by 


velocity  (VQ)  is  standardly  (Reference  5  ) 


(A. 10) 


where  G  is  the  ventiliation  coefficient  (G  =  .711)  and  p  is  the 
trench  gas  density. 

The  mass  outflow  rate  per  unit  length  (mQ)  is 

4Gx£cos20  h 

mo  =  PVow  =  - - -  [2(P-P0)p]  (A.  11) 

If  the  venting  velocity  exceeds  the  local  trench  sound 
speed  (c)  ,  the  venting  flow  is  assumed  to  choke  at  speed  c. 

In  this  case 

m0  =  pew  =  2pcxrcos0,  V0>c  (A. 12) 


The  radius  ( R)  of  an  expanding  cylindrical  tunnel  was 
modeled  by  letting 

2  -  X_R 
W 


r2  -  Ro~  Cl+-p) 


(A. 13) 


The  change  in  total  gas  density  in  the  trench  due  to 
venting  and  wall  expansion  can  then  be  expressed  by 


P  =  P'vent  +  ^ 


expansion 


where 


and 


p'' 


m. 


vent 


ttRA 


pi  =  .  iei 

expansion  ^ 


(A. 14) 

(A. 15) 

(A. 16) 
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